Numerical simulation of moving contact line problems using a volume-of- uid method

Moving contact lines are implemented in a volume-ofuid scheme with piecewise linear interface construction. Interfacial tension is treated as a continuous body force, computed from numerical derivatives of a smoothed volume-ofuid function. Two methods for implementing the contact angle condition are investigated. The rst extrapolates the volume-ofuid function beyond the ow domain, on the basis of the condition that its gradient is perpendicular to the interface and the normal to the interface at the wall is determined by the contact angle. The second method treats the problem as a three-phase situation and mimics the classical argument of Young. It is found that the latter approach introduces an arti cial localized ow, and the extrapolation method is preferable. Slip is a crucial factor in the spreading of contact lines; the numerical method introduces slip at the discrete level, e ectively introducing a slip length on the order of the mesh size.